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Abstract : Main effect plans orthogonal through the block factor (POTB) have been defined and 
a few series of them have been constructed in Bagchi (2010). These plans are very closely related to 
the ‘mutually orthogonal balanced nested row-column designs’ of Morgan and Uddin (1996) and many 
other combinatorial designs in the literature with different names like ‘BIBDs for two sets of treatment’, 
‘Graeco-Latin designs’ and ‘PERGOLAs’. In fact all of them may be viewed as POTBs satisfying one or 
more additional conditions, making them ‘optimal’. However, the PERGOLAs are defined to satisfy an 
additional property, without which also it is optimal. Interestingly, this additional property is satisfied by 
all the hitherto known examples of POTBs, even when their definitions do not demand it. 

In this paper we present direct and recursive constructions of POTBs. In the process we have con¬ 
structed one design which seems to be the first example of an ‘optimal’ two-factor POTB which is not a 
PERGOLA (see Theorem 13.II) . 

1 Introduction. 

Preece (1966) constructed ‘BIBDs for two sets of treatments’. Subsequently several authors constructed 
similar combinatorial objects. Among these, the ones relevant to the present paper are ‘balanced Graeco- 
Latin block designs’ of Seberry (1979), ‘Graeco-Latin designs of type 1’ of Street (1981) and ‘Perfect 
Graeco-Latin balanced incomplete block designs (PERGOLAs)’ of Rees and Preece (1999). 

Morgan and Uddin (1996) considered main effects plans (MEPs) on a nested row-column set up and 
proved the optimality of ‘mutually orthogonal balanced nested row-column designs’. They also discussed 
the constructional aspects of such designs. Unfortunately, the relevance of these results in the context of 
blocked MEPs was overlooked by later authors studying blocked MEPs like Mukerjee, Dey and Chatterjee 
(2001) and Bagchi (2010). Optimal Main effects plans for three or more factors on non- orthogonal blocks 
of small size were obtained in Mukerjee, Dey and Chatterjee (2001). Bagchi (2010) defined and studied 
main effect plans orthogonal ‘through the block factor’ (POTB). 

In the present paper we first note the relation between POTBs with optimality property (termed 
balanced POTB) and the combinatorial objects considered by earlier authors mentioned above. Next we 
construct a few series of POTBs. We also present a recursive construction by which the number of factors 
is multiplied, keeping the block size unchanged, thus yielding multi-factor POTBs from PERGOLAs and 
other similar two-factor designs. 

We note that a PERGOLA is a balanced POTB with an additional condition. It is interesting that all 
the balanced POTBs available in the literature (with different names) do satisfy this condition. (Table 1 
of Rees and Preece (1999) shows that such designs are plentiful). One would, therefore, suspect that this 
condition is implicit in the definition. We have, however, found a balanced POTB which does not satisfy 
this condition. [See Theorem 13.1| . 
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2 Preliminaries 


In this section we present the definition of a balanced POTB. We also list related combinatorial objects 
existing in the literature with various names. 

Definition 2.1 By a block design with v treatments and b blocks of size k each we mean an incidence 
structure represented by a v x b matrix N having constantcolumn sum k. 

With any such block design D one associates the graph G(D) with the treatments of D as vertices, two 
treatments being adjacent in G(D) if there is a block containing both the treatments. One says that D is 

connected if the graph G(D) is connected in the usual sense. 

V(b, k, v ) will denote the class of all connected block designs with v treatments on b blocks 
of size k each. 

Definition 2.2 Consider an (m + 1) x n array A such that the entries of the ith row are elements of a 
set Si of size Si, i = 0,1, ■ ■ ■ m. This is said to be a main effect plan (MEP) for m+ 1 factors, say, 
Fo,F\, ■ ■ ■ F m on n runs. The ith row corresponds to the factor Fi and we say that Ft has Si levels. 

For 0 < i, j < to, let M,.y be the Si x Sj matrix such that the rows and columns of M,^ are indexed 
by S i and Sj respectively and the (p, g)th entry of My is the number of columns of A in which 
the ith and jth entries are p and q respectively, p € Si, q £ Sj. My is said to be the Fi versus 
Fj incidence matrix. 

Now, suppose n = bk, where b and k are integers. A blocked MEP (laid out on b blocks 
of size k each) is an (m +1) x bk array in which So is the set of integers {1, 2 ■ ■ • 6} and each integer 
j £ So appears exactly k times in the last row of A. The O-th row is said to correspond to the 
“block factor” , which is represented by B (and not Fq). In this case, the incidence matrix of the 
ith row versus the last row is denoted by Mis,l < i < m. A blocked MEP is said to be 
‘connected’ (borrowing a term from the theory of block designs) if each M,s is the incidence matrix 
of a connected block design. It is said to be symmetric if si = s% ■ • • = s m . A symmetric MEP with 
Si = s for every i is also referred to as an MEP for an s m experiment. 

In the applications, there are n experimental units, which are classified into homogeneous classes or 
blocks. These units are used to study the effects of m factors, the ith one having Si ‘levels’. Typically, an 
experimental unit, say in the jth block, receives a ‘level combination’ say x = (xi, ■ • ■ x m ), i.e. the level Xi 
of Fi is applied on that unit, i = 1,2, ■ • - to. This information is stored in the column vector ( j, x\, ■ ■ ■ x m )'. 
The array A consists of all such column vectors. 

Definition 2.3 [Bagchi(2010)J The ith and jth factors of a blocked MEP p are said to be orthogonal 
through the block factor (OTB) if 


MjB(MjB)' = kM.ij. (2.1) 

If every pair of factors of a plan p is orthogonal to each other through the block factor, then p is said 
to be a plan orthogonal through the block factor (POTB ). 

Remark 2.1: What is the utility of condition (1 2.1 I) ? This condition guaranties that for the inference 
on a factor Fi of a POTB one has to look at only its incidence with the block factor (i.e. Mib) and forget all 
other treatment factors. Thus, the performance of a POTB p regarding the inference on the ith treatment 
factor depends only on the incidence matrix Mis -We present a more precise statement in the next theorem. 
We omit the proof which can be obtained by going along the same lines as in the proofs of Lemma 1 and 
Theorem 1 of Mukerjee, Dey and Chatterjee (2001). [See Shah and Sinha (1989) for definitions, results 
and other details about optimality]. 
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Theorem 2.1 Suppose a connected POTB p* satisfies the following condition. For some non-increasing 
optimality criterion </>, MiB is the incidence matrix of a block design d which is (f> -optimal in the class 
of all connected block designs with Si treatments and b blocks of size k each. Then, p* is (f>-optimal in the 
class of all connected m-factor MEP in the same set-up for the inference on the ith factor.. 

In particular, using the well-known optimality property of a BIBD we get the following result. 

Corollary 2.1 Suppose p* is a connected POTB. Suppose further that Mm is the incidence matrix of a 
BIBD. Then, p* is universally optimal in the class of all m-factor connected MEP in the same set up, for 
the inference on the ith factor. 

In view of the above result, we introduce the following term. 

Definition 2.4 A connected POTB is said to be balanced if e ach of its factors form a BIBD with 
the block factor, that is Mm is the incidence matrix of a BIBD for each factor Fi. 

We now present a small example of a balanced POTB. This has two factors, each with four levels 
0 , 1 , 2 ,3 on six blocks of size two each. 

Blocks \ B\ | f ?2 | \ B 4 B 5 ^6 

F\ j 0 T 3"1 6 3"1 I —j 0 I 1 3 2~ 

F 2 |13i02i21|30|32|01- 


Next we list a few combinatorial designs and note their relation with balanced POTBs. 

(a) Balanced Graco-Latin block design defined and constructed in Seberry (1979) are balanced 
POTBs with two factors. 

(b) Graco-Latin block design of type 1 of Street (1981) are also two-factor balanced POTBs having 

M 12 = J. 

(c) Perfect Graeco-Latin balanced incomplete block designs (PERGOLAs) defined and dis¬ 
cussed extensively in Rees and Preece (1999) are two-factor balanced POTBs having 

(i) si = s 2 = s, say and 


M 12 M' 12 = M' 12 M 12 = fl s + gJ s , f, g are integers. ( 2 . 2 ) 

(d) Mutually orthogonal BIBDs defined and constructed by Morgan and Uddin (1996) are multi¬ 
factor balanced POTBs. 

Here I n is the identity matrix and J n is the all-one matrix of order n. 

Remark 2.2: The definition of neither balanced Graco-Latin block designs nor of mutually orthogonal 
BIBDs include condition (1 2.2 I) . However, it is interesting to note that all these designs constructed so 
far do satisfy this condition. [See theorem 13.5] . 

3 Constructions for symmetric POTBs 

Now we present a few constructions of POTB’s. Each of these constructions is in terms of some group G 
of order g (which is the additive group of the field V in Theorem 13.41) . 

A block will consist of k plots or runs represented by columns. By adding an element 
u £ G to a block we mean adding u to the level of every factor in every run of the block. By 
developing an initial block we mean generating g blocks by adding distinct elements of G to 
the initial block. 

Let N denote the set of integers modulo n and N + denote N U { 00 }. 
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Theorem 3.1 Let n be a positive integer > 5. 

(a) Then there exists a POTB with three factors Fq, F 4 , F 2 each having n + 1 levels on b = 6n blocks 
of size two. 

(b) In the case n = 5, we get a Balanced POTB, which is not a PERGOLA. 

Proof : (a) Let N + be the set of levels for each factor. The initial blocks Bij, i = 1, 2, j = 0,1, 2 are 
as follows. Here addition in the suffix of F is modulo 3. 



Block 

Bij 

Block 

B2.i | 

F 0+j 

oo 

0 

oo 

o 

Fi+j 

1 0 

1 

1 0 

2 

F 2 +j 

1 

1 

1 1 

2 


That the design satisfies the required property follows by straightforward verification, 
(b) Let n = 5. One can verify that the incidence matrices satisfy the following. 


Mij 


M iB (M iB y 


0 2 2 2 2 2 

2 2 2 1 1 2 

2 2 2 2 1 1 

2 1 2 2 2 1 

2 112 2 2 

2 2 112 2 

8 Iq + 2 Jq, i = 1,2, 3. 


i,j = 0,1, 2 and 


(3.3) 


(3.4) 


We see that each M iB is the incidence matrix of a BIBD. Thus, by definition ^. 4l it is a balanced POTB. 
But neither of M i3 s satisfy (1 2.2 ll . as is clear from (1 3.3 I) . Thus, the two-factor balanced POTB obtained 
by ignoring any one of the factors is not a PERGOLA. □ 


Theorem 3.2 Suppose n is an integer > 5. Then there exists 

(a) a POTB with two n— level factors on 2 n blocks and 

(b) a POTB with four n—level factors on An blocks of size 2 each. 

(c) We get a balanced POTB in the case n = 5 in series (a). Further, in the case n = 10 in series (b) 
we get a POTB which is E-optimal for the inference on each factor. 


Proof : The set of levels of each factor is N. Let a,b £ N. 

(a) We present initial blocks B 4 and B 2 below. 

Block B\ | Block B 2 \ 

Fi | a —a | b — b \ . 

F 2 | b — b | —a a 

(b) We present the initial blocks Bi , Z = 1, - - • 4 below. 


Blocks 

Bi 

1 

b 2 



B 3 


1 B 4 


Fi 

0 

a 

a 

—a 


0 

b 

1 -b 

b 

f 2 | 

a 

—a 

0 

—a 


-b 

b 

o 

b 

f 3 | 

0 

b 

b 

-b 


—a 

0 

a 

—a 

f 4 : 

b 

-b ; 

0 

-b 


a 

—a 

a 

0 


That these initial blocks generate POTBs can be verified by straightforward computation. 

(c) For n = 5, taking a = 1, b = 2 we get a balanced POTB. 

For n = 10, we take a = 1 and 6 = 3. Then for every i = 1 , • • • 4, MiB is the incidence matrix of a group 
divisible design with five groups, the jth group being the pair of levels {j, j + 5} j = 0, • • - 4, satisfying 
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Ao = 0 and Ai = 1. This plan is, therefore, E-optimal for the inference on all the four factors by Takeuchi 
(1961). □ 


Theorem 3.3 (a) There exists a symmetric POTB with four n—level factors on 4 n blocks of size 2 each, 
whenever n > 9. We get a balanced POTB in the case n = 9. 

(b) There exists a symmetric POTB with four factors each having n + 1 levels on 6 n blocks of size 2 
each, whenever n > 7. 


Proof : (a) The set of levels for each factor 
are as follows. 

Blocks Bi B 2 \ B 3 

is N 

. Let a, b,c,d € N. The initial blocks Bi , l = 1, • • • 4 

I Ba 

Fi 

1 

a 

—a b 

-b 

c 

—c 

| ~d d 

f 2 

1 

b 

—b —a 

a 

-d 

d 

—c c . 

f 3 

1 

c 

—c d 

-d | 

—a 

a 

b -b 

Fa 

1 

d 

—d —c 

c 

b 

-b 

a —a 

One can 

easily verify that these initial blocks generate a symmetric POTB with the given parameters. 

By taking a 

= 

1 ,b = 

2, c = 3 and d 

= 4 in 

the case n 

= 9, we get a balanced POTB. 

(b) The set of levels for each factor is 

N+. 

Let 

x,b,c € N. The initial blocks Bi, l = 1, • • • 6 are as 

follows. 








Blocks 

1 

Bi 

1 b 2 

1 

b 3 


| Ba \ B 3 \ Bq 

Pi 

1 

0 

oo | a 

—a 

b 

~b 

\ c —c | a —a | a —a 

f 2 

1 

a 

-a | 0 

oo 

c 

—c 

—b b | a —a \ —a a 

f 3 

1 

b 

—6 j c 

—c 

0 

oo 

a —a | —c c | — c c 

Fa 

1 

c 

—c b 

-b 

a 

—a 

0 oo | — c c | c — c 


That the design satisfies the required property follows by straightforward verification. □ 

Next we construct a series of balanced POTBs using finite fields. We first introduce the following 
notation. 

Notation 3.1 (i) Let |J denote an union counting multiplicity. 

(ii) For a set A and an integer n, let nA denotes a multiset in which every member of A occurs n 
times. 

(in) For subsets A and B of a group ( G , +), 

A — B = {a — b : a £ A, b € B}. 

Notation 3.2 (i) v denotes an odd prime or a prime power, written as v = mf + 1. V denotes the Galois 
field of order v. Further, V* = V \ {0} and V + = V U {oo}. 

(ii) a denotes a primitive element ofV. 

(Hi) fd = a m is a generator of the subgroup Co of order f of (V *,.). 

(iv) Cq, Ci, ■ ■ ■ C m -i are the cosets of Cq in (V*,.). 

(v) ( i,j ) = the number of ordered pairs of integers (s,t) such that the following equation is satisfied in 
V. [ This notation is borrowed from the theory of cyclotomy] 

1 + a s = a*, s = i,t = j (mod to). 

We need the following well-known result. [See Hall (1986), for instance]. 
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Lemma 3.1 Suppose m = 2. Then the following hold. 

(a) —IS Cq (respectively C\) if f is even (respectively odd). 

(b) If f is even, then a - 1 S Cj => a -1 - 1 £ Ci+i ■ 

(c) If f is odd, then a - 1 £ Cj 1 - a -1 £ Cj+i- 
Here + in the suffix is modulo 2. 

We present the following well-known results for ready reference. [See equations (11.6.30), (11.6.40) and 
(11.6.43) of Hall (1986)]. 

Lemma 3.2 The differences between the cosets Cj’s of V* can be expressed in terms of the cyclotomy 
numbers (p,q)’s as follows. 

f m— 1 

U (k~j,i~j)Ck tfj^i 

k =0 

m—1 

/{ 0 }U U [k-mCk ifj = i 

k—0 

The following cyclotomy numbers are known for the case to = 2. 

Case 1: f odd. (0,0) = (1,1) = (1,0) = (f-l)/2, (0,1) = (f+l)/2. 

Case 2: f even. (0,0) = f/2 -1, (0,1) = (1,0) = (1,1) = f/2. 

A series of two-factor balanced POTBs : 


Theorem 3.4 Suppose v is an odd prime or a prime power. Then there exists a balanced POTB for a 
(v + l) 2 experiment on b = 2 v blocks of size (v + l)/2 . 


Proof : We write v = 2/ + 1. The set of levels of each factor is V U {oo}. The plan is obtained by 
developing the following initial blocks Bq and B\ presented below. 

Case 1 : f is even. 


Bo 


oo 1 (3 ••• pf- 1 
0 a a/3 ■ ■ ■ afif- 1 


and B\ 


0 1 p 

oo a -1 a -1 /3 


Case 2 : f is odd. 

Block Bo is as in case 1, while Block Bi is as follows. 


Pf- 1 1 

a-'pf - 1 


0 a 1 a 1 /3 • • • a 1 /3 f 1 
oo 1 P ■ ■ ■ pf- 1 


Clearly block size is / + 1 = (v + l)/2. To show that the plan satisfies the required property, we have 
to show that 

(a) the plan is a POTB and (b) each factor forms a BIBD with the block factor. 

Condition (b) follows from Lemma I5~21 So, we prove (a) . 

let us use the following simplified notation M = ((to^)) for M 12 and A for M 1 b(M 2 b) / . We note 
that niij is the total number of plots (runs) in which the level combination (i,j) appears, while aij is the 
number of blocks in which F\ is at level i and F% at level j, (in same or different plots). 

We shall show that 


M = J-I and (3.5) 

A = (/ + 1)(J-I) (3.6) 

We begin with M. It is clear that mu = 0,*£ V + and m a = mi t00 = 1 ,i £ V. 
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We, therefore, assume i ^ j, i,j £ V. Let u = j — i. Then, is the number of times u appears in 
the multiset 


f (a — l)Co |J( a 1 — l)Co if f is even 
\ {a — 1)C 0 1 |(1 — a~ 1 )Co iff is odd 

The relations above imply (1 3.5 II in view of Lemma 13.11 

Now we look at A. Clearly, A ti = 0, i £ V + . Further, for every i £ V, moo,* is the replication number 
of i in the block design generated by the initial block {0} U C\. Similarly, mi i00 is the replication number 
of i in the block design generated by the initial block {0} U Co if / is odd and {0} U C\ otherwise. Thus, 

nToc,* — o — f A 1, z £ V. 

Next we consider i ^ j, i,j £ V. Let u = j — i. Then, aij is the number of times u appears in the 
multiset 

of (({0}UC 1 )-C o )U(Ci-({0}UC'o)) iff is even 

\ (({0}UC 1 )-Co)U(Co-({0}UC 1 )) iff is odd. 

These, together with Lemma [3.21 and (a) of Lemma 13. II imply the equation next to (1 3.5 I) . □ 

Now we present the series of balanced POTBs available in the literature, together with two newly 
constructed balanced POTBs. 

Notation 3.3 v denotes an odd prime or a prime power of the form v = mf + 1. 


Table 3.1 : Balanced POTBs 


No. 

The expt. 

m 

ff of Blocks 

Block size 

Inc. matrix (My) 

Reference 

1. 

(d + 1)x» 

2 

2 v 

/ = («- l)/2 

J 

Seberry (1979) 

2. 

{V + l) 2 

2 

2 v 

/ + 1 = (u + l)/2 

J 

Street (1981) 

3(a). 

v m 

m 

mv 

/ 

C J-i) 

Morgan and Uddin(1996) 

3(b). 

V* 

tg 

mv 

hf, h<g 

h(J — I) 

Morgan and Uddin(1996) 

4. 

v f 

m 

mv 

1 + hf, h < m 

(?n — h)I + hJ 

Morgan and Uddin(1996) 

5. 

(V + l) 2 

2 

2 v 

/ + ! 

J-I 

Theorem 13.41 

6. 

(5 + l) 3 

- 

30 

2 

riven in(| 3.3 |) 

Theorem 13.11 

7. 

9 4 

- 

36 

2 

J-I 

Theorem 13.31 


In 3(b) above t is a factor of m,g — m/t and h is an integer < g. 


Using the information in Table 3.1, one may verify the following result. 

Theorem 3.5 Every two-factor balanced POTB obtained from an existing multi-factor balanced POTB , 
except the one with n = 5, constructed in Theorem \3.1[ is a PERGOLA. 

One may also look at Table 1 of Rees and Preece (1999) for many more examples of PERGOLAS. 

A recursive construction 

Notation 3.4 An orthogonal array with m rows, n columns, k symbols and strength 2, will be denoted by 
OA{n, m, k, 2). 

Theorem 3.6 Suppose there exists a balanced POTB with f factors on b blocks of size k each, with f < k. 
If further an OA(k 2 , m + 1, k, 2) exists, then a balanced POTB with mf factors on bk blocks of size k each 
also exists. 
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The proof of this theorem is based on the following lemma. 

Lemma 3.3 Consider a set of k runs of a plan for an experiment with /(< k) factors, such that no level 
of any factor is repeated. If an OA(k 2 , m + 1, k, 2) exists, then there exists an MEP with f classes of m 
k—level factors on k blocks of size k each with the following property. Every factor is orthogonal (w.r.t. 
the block factor) to every factor of a different class. 

Proof : Let D denote the given set of runs. Let F = {P, Q, ■ ■ ■ } denote the set of f factors of D. For 
each P £ F, let Kp denote the set of levels of P appearing in D. Let K = {1, • - - fe} denote the set of 
symbols of the given OA (O), say. For every P £ F, let Lp denote the following one-one function from K 
to Kp. [By assumption, size of Kp is k for each P\. 

Lp(i ) = j, i £ K,j £ Kp, if P has level j in the ith run of D. (3-7) 

Let us arrange the columns of the given OA (O) as 

(0)=[A 1 

such that the 1st row of Aj consists of the symbol i repeated k times, 1 < i < k. Let Ai denote the 
m x k array obtained from Aj by deleting the 1st row. Thus, every member of K appear exactly once in 
every row of each Ai, i = 1, • • • k. 

We now construct D*, the reqd MEP. For each factor P of D, there will be m factors Pi, • • • P m in D*, 
each of which will have Kp as the set of levels. 

For i £ K , let us fix Ai and a factor, say P of D. If the j th column of Ai is (si, • • • s m ), s u £ K, then 
in the j-th plot of the Pth block of D*, the factor P t will have level Lp(st), t = 1, 2, • • • m, where Lp is as 
in (1 3.7 I) . Doing the same for all the factors and varying j over {1, 2, • • • k} we get a block of D*. Finally 
varying i over K we we generate the k blocks of the reqd MEP. 

We now show that the MEP D* satisfies the required property. We fix two factors, say Pi and Qj, i A j 
and an ordered pair of levels, say (u,v),u £ Kp, v £ Kq. From the construction the following is clear. 
In every block there is a plot in which P* is at level u and a plot where Qj is at level v. Moreover, there 
is exactly one block in which these factors are set at these levels in the same plot. Thus, the factors Pi 
and Qj are mutually orthogonal through the block factor. We see that if P = Q, then also the argument 
above holds. Thus Pi is orthogonal to Pj, j A i. However, Pi may not be orthogonal to Qi. We, therefore 
form the classes as Cf = {Pi,Qi, ■ ■ ■ }, P,Q £ F, i = 1, 2, • • • m. Now the factors satisfy the orthogonality 
condition of the hypothesis.□ 

Proof of the theorem : Let D denote the given POTB. Let Ai,i £ K be as in Lemma. For every 
block of D we construct a an MEP following the method described in the proof of the lemma above. Let 
the resultant MEP be named D*. By Lemma T3.31 every pair of factors belonging to different classes are 
orthogonal w.r.t. the block factor. Further, since the pair of factors P, Q are mutually orthogonal w.r.t. 
the block factor in D, it follows that for every i £ K, the factors Pi and Qi are also mutually orthogonal 
w.r.t. the block factor in D*. □ 

Remark 3.1: If we look at the restriction of D* to one factor, say Pi, we see that it is nothing but k 
times repetition of each block of the block design obtained from the restriction of D to the factor P. 
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